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£P! Abstract 
P 

In this paper we derive the finite size corrections to the energy eigenvalues of the energy for 
2D dimers on a square lattice. These finite size corrections, as in the case of Critical Dense 
Polymers, are proportional to the eigenvalues of the Local Integrals of Motion of Bazhanov 
Lukyanov and Zamolodchikov for central charge c = —2. This sheds more light on the status 
of the Dimer model as a conformal field theory with this value of the certral charge. 

1 Introduction 

In [1] Lieb proposed a transfer matrix for the Dimer model, recently in [2] Rasmussen and 
Ruelle exploited this transfer matrix description to give a partitioning of the dimer config- 
urations into sectors which are closed under the action of the transfer matrix. In the case 
of open spatial boundary conditions each sector was found to be associated with a repre- 
sentation of the Virasoro algebra, since the partition function of a sector yields a Virasoro 
character in the continuum scaling limit. The characters can be organized consistently with 
a CFT of central charge c = —2. 

The selection rules for the eigenvalues of the transfer matrix are reminiscent of those for 
Critical dense Polymers [3], [I], since the eigenstates of the transfer matrix can be labelled 
by two column diagrams as well. 

In [5] Ruelle et al. analyzed the finite size corrections to the energy levels for critical dense 
polymers, and found that the first correction can be described by conformal perturbation 
theory made with operators belonging to the tower of the identity. The analysis of finite size 
corrections, however had already been carried out to all orders in [B] albeit not by means 
of conformal perturbation theory. In this work it was found that each finite size correction 
comes with a coefficient which is proportional to the eigenvalues of the BLZ |7] Local Inte- 
grals of Motion (IOM), from which it follows immediately that the ratio of the amplitudes 
for finite size corrections is universal. 

In this work we shall compute all the finite size corrections for the energy levels of the Dimer 
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model, and show that the same result of [6J works in this case as well, i.e. that the finite 
size corrections are proportional to the eigenvalues Local IOM of [7], as computed in [6]. 
We believe that this result is interesting since it yields universal quantities of c = — 2 CFT 
without using the value of the central charge as an ansatz which is checked then by con- 
sistency of the results. This work is organized as follows: in sections 2,3,5 we very briefly 
review the results of [2J, whereas in section 4 we derive the exact form of the finite size 
corrections for Dimers. 



2 Dimers and Transfer Matrix 

In this section we shall briefly review the results of PQ[2] concerning the definition of the 
Dimers model and of a suitable transfer matrix. 

A Dimer model on a square lattice of size N x M is a model where the lattice is covered 
by elementary tiles covering exactly 2 sites, the Dimers being not allowed to overlap. Each 
Dimer can be either orizzontal or vertical, to each dimer we associate a weight which can 
be rescaled so that a vertical Dimer has weight 1 and a horizontal dimer has weight a. The 
partition function of such a model can be written as 

Zn,m= aH 

configurations 

being h the number of horizontal Dimers. 

There is a mapping between Dimers on a square lattice and configurations of up and down 

arrows on the same lattice: and up arrow at site i means that there is a vertical Dimer 

covering that site and its northern neighbour, a down arrow represents the absence of such 

a dimer at site i. A row of N sites can have 2 N Dimer-arrow configurations. 

It is then intuitive to think of each arrow as an element of a canonical basis (1, 0), (0, 1) of 

C 2 , a row configuration then lives in the iV-fold tensor product of C 2 with itself. 

We now want to build a Trasfer Matrix in terms of action of tensor products of pauli matrices, 

we thus define 

<7t = l®l®...<8><7®...®l (2.2) 
where a is in position i in the tensor product and is one of the following matrices: 

-(:;)■ ^-(:s). **-(?;) <-> 

The Transfer Matrix is then defined as the product 

T = (2.4) 

where 



N 



(2.5) 



8=1 
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and 

N-l 

V 3 = H(l + aa-ar +l ) (2.6) 
i=i 

The action of V\ flips all the arrows in a given row, whereas in the action of V3 each factor 
either lives the local arrow unchanged by acting with the term 1 or produces a horizontal 
dimer by reversing a pair of neighbouring up arrows by acting with anc ^ yielding 

a weight a contribution. The product V3V1 generates all possible ways of going from a 
prescribed set of arrow configurations by creating a horizontal dimer whenever possible and 
correctly handling incoming vertical dimers. 

The partition function on a strip with periodic boundary conditions in the M direction is 
then 

Z = TrT A/ (2.7) 

It happens that such a partition function vanishes if both N, M are odd, therefore we shall 
assume M to be even, this also enables us to use in our analisys T 2 instead of T itself. 
The transfer matrix T 2 can be diagonalized by a Jordan- Wigner transformation followed by 
a Fourier transform, and its eigenvalues are of the form, for N = 2L even: 

A = jg ^ 1 + a* + „ B in (*||^)) 2<1 " <i "" i) (2.8) 

whereas for odd N = 2L + 1 

A = rj(^l + c? sin* + asin (^)) ' ' (2.9) 

where €j, \ij independently take the values 0, 1 to such a pair (e 3 -| it is possible to associate 
a 2-column diagram V = (ji, . . . ,ji\ii, . . . , i r ) there the js are those for which ej = 1 and 
similarly the is are those for which /ij = 1. 

The 2-column diagrams T> will satisfy certain selection rules which we will describe in the 
following. 

The energy eigenvalues are obtained from the As by taking the logarithm 

E(V) = - l -\og\ (2.10) 

3 Sectors and Variation Index 

In this section we show how to partition the space of states of the Dimer model into sectors 
labelled by the eigenvalues of an operator called variation index [2]. 
Following [2] let us introduce the operator 

1 - 

i=i 
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This operator is digonal in the arrow basis and has eigenvalues 



+ (3.2) 

we notice that the eigenvalues are integers for N even, and half integers for N odd. 

Let E v be the vector space spanned by states having a fixed v, it can be shown that the 

generating function for the dimensions dimE v is: 

JV 

dimE v t v = (y/i + -j^ (3.3) 

from which it follows that: 

dim^ = dim£_„ = \^ N N J (3.4) 

we now want to study how the transfer matrix T = V3V1 acts on the spaces E v , it is 
not difficult to see that V anticommutes with V\ and commutes with V3, it follows that it 
anticommutes with T, thus mapping 

T : E v ^ E_ v (3.5) 

it follows that T 2 commutes with V and that each linear space E v is invariant under the 
action of T 2 , whereas the T- invariant subspaces are of the form: 

E v © E_ v (3.6) 

we call E v © £L„ (and E$) T-sectors and E v T 2 -sectors. 

4 Finite Size Corrections 

In this section we derive the finite size corrections for the Dimer model and show that they 
are proportional to the eigenvalues of the integrals of motion of [7j. 

4.1 Even N 

Let the size of the system be N = 2L, one then has that the logarithm of eigenvalues of the 
transfer matrix is given by 

.7r(2fc-l)v\ A / . ,n(2k-l) 



E = ^arcsmh(asm j - ^ axcsmh [a sin [ —-— ) j (4.1) 

k == 1 

if we now consider the following expansion 



00 



we have that 

E = Y(-l) k , ^* , a 2fc+1 (Y-Y) sin 2fe+1 ( < 2m - l } \ 
^ y ' (2k + l)k\ \^ V 22L + 1 ; 

k=0 v ' m£T> m=l v 7 



(4.3) 



we now have that 

y.. at+ i M2m-l) N 2L+l (l) t 1 ~ / T s a + i fl 2 , +2 (l) 

E 8m I W —(15 " 2 " £ (amJ pTT2)! C « (44) 



oo 22+1 

^ 2n+1 W = E c ^(|TI)! (45) 

being 

(2(fc - n) - If 1+1 (4.6) 



^^£(-<r> 

i — n \ / 



k=0 

it is then straightforward to obtain that 

V" _ V" ^ Q ir,2*+i / 7[(2m-l) v _ 2L + 1 (l) fc 1 
2^ 2^j sm ^ 2 (2L + l) j " tt (§) fc + 2 + 



+ 2^{ 2 L + l) (2l+l)\2U+A^ Zm ij + 21 + 2 

(4.7) 

we now notice that 

. o22+lc> m. 9/-I 

EP»-ir + ^ (4 - 8) 

where the /2Z+1 are the eigenvalues of the integrals of motion of BLZ [7] in the sectors of 
conformal weights A = — §, §> 

Now, summing up the series for the constant and divergent part gives the bulk and boundary 
free energy contribution, and collecting a fixed order of 2L + 1 in the sum gives the following 
expansion from which we read off all the finite size corrections: 

2L + 1 ^ / 7r P,(a) 

E=— /bull ( Q ) + /bou ( a )+£(_) (2i+1) ,' 2 U (t+1) WP) (4.9) 

where the polynomials Pi (a) are defined as follows 

^=h-^jS^ c '* a2M (4 - io) 
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the first few of these polynomials are given by 

P (a)=a (4.11) 

P L (a) = -a(l + a 2 ) (4.12) 
P 2 (a) = a(l + 10a 2 + 9a 4 ) (4.13) 

where 

7T 

/buik(«) = — / dt arcsinh(a sin t) (4-14) 
Jo 

and 

/bou(a) = -arcsinha (4.15) 

it can be checked that the expansion of /bulk (a) in powers of a is correctly reproduced by 
the coefficients of the expansion for the divergent part obtained here. 



4.2 Odd N 

Let the size of the system be N — 2L + 1, one then has that the logarithm of eigenvalues of 
the transfer matrix is given by 

E = arcsinh (a sin ( ^ + ^ ) ) - ^ arcsinh (a sin ( ^ + ^ ) ) (4.16) 

we now have that 

^ Sm l(2tT2) J " — (|U " 2 " 2- 12LT2 J (27T2)! Cw (4 ' 17) 
it is then straightforward to obtain that 



(£-£H' +1 (^) 



urn . _ 2L + 2(l) fc 1 
(2TT2) j_ ^ + 2 + 



meD m=l v y v 2 



+ E ( 2L + 2) (2/ + l)Wi ( E( 2m ) 2m + 2/ + 2^ } ) 

i=k \ ' j meV 

(4.18) 

we now notice that 



where the / 2 ;+i are the eigenvalues of the integrals of motion of BLZ in the sectors of 
conformal weights A = 0, 1, 3, 6, . . .. 
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Now, summing up the series for the constant and divergent part gives the bulk and boundary 
free energy contribution, and collecting a fixed order of 2L + 2 in the sum gives the following 
expansion from which we read off all the finite size corrections: 

E = 2 ±±l Ma) + /bou(Q ) + g (^)"' (i , + ff (itl) Mi)) (4.20) 



5 Characters and Selection Rules 

In this section we review the selection rules for the two-column diagrams T> labeling the 
eigenstates of the integrals of motion appearing in the finite size corrections. 
This results follow the exposition given in [2J, and the conventions for the irreducible fini- 
tized characters ch^P (g) of introduced in the work on critical dense polymers [3] . 



5.1 Even N 

We first of all notice that the variation index v can be conveniently parametrized as: 



v = £(^ 



(5.1) 



we then recall that by neglecting the divergent and constant part the energy levels behave 
as: 



E 



an 



c 1 



£a - b 



iV + l v 24 8 ' ^ w T 



the finitized partition function of the sector E v then reads 



being 



v&u. 



Mar 

q = e N+i 



(L) 



(5.2) 



(5.3) 



(5.4) 



where U% is the set of all two column configurations of variation index v and maximal height 
L. 

It is possible to see that the finitized character of the sector E v is also expressible as: 



where 



denotes the q-Binomial coefficient and 



2L 
L — v 



(5.5) 



(2r - s) 2 - 1 



(5.6) 



is the conformal weight in the Kac table. 

Now, by introducing the finitized irreducible characters 



ch S +1) (?) = Q 24 



+A r , 2 



l-q 



2r 



l-q 



N+2 



N 



N + 1 
r + 1 



one can write 



LorL+l 



and in the limit N — > oo: 



zi N) = E ch S +1) (^) 

r=|u|+l,6j/2 
oo 

r=\v\+l,by2 



(5.7) 



(5.8) 



(5.9) 



where ch r tS (q) are the irreducible characters of conformal weight A r jS in the continuum limit. 
The above expression can be simplified to: 



z v {q) = 



77(g) 



(5.10) 



where i](q) is the Dedekind rj function. 



5.2 Odd N 

Let us introduce the excess parameter w defined as: 

L 

i=i 

since the variation index v is in this case a half integer its values are related to the excess 
parameter by 

w = v±^ (5.12) 

which thus has values 

w G {—L, —L + 1, . . . , L — 1, L} (5.13) 
again by neglecting the divergent and constant part the energy levels behave as: 

E ( h> ?) (5.14) 

jev 

the finitized partition function of the sector E v then reads 

Z^(q) = q-^ E q L * vi (5.15) 
v&u {L \ uu (L \ 



8 



being 

Mock 

q = e N + 1 (5.16) 

where U% is the set of all two column configurations of excess parameter w and maximal 
height L. 

It is possible to see that the finitized character of the sector E v is also expressible as: 



2L + 1 



Now, by introducing the finitized irreducible characters 



ch^ +1) ( ? ) =q 24- 



-Ar.l W 



iV + 1 



_ r+ 1 

L 2 ' ~ 2 -lg 



one can write 



(5.17) 



(5.18) 



(5.19) 



and in the limit iV — > oo: 

oo 

r=|«|+i 

The above expression can again be simplified to: 



(5.20) 



77(g) 
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